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Nomenclature 

fraetion  of  molecules  reflected  diffusely  from  a  surface 

area,  m^ 

diameter,  m 

ratio  of  specific  heats 

specific  impulse,  sec 

Boltzmann’s  constant  =  1.3806503  x  10'^^  J/K 

nominal  length  of  the  diverging  section  of  the  nozzle,  m 

nominal  length  of  the  diverging  section  wall  of  the  nozzle,  m 

molecular  mass,  kg 

viscosity,  kg/m  s 

Mach  number 

unit  vector  normal  to  a  surface 
density,  kg/m^ 
pressure.  Pa 
Reynolds  number 
orifice  thickness,  m 
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T 

3 

T 

e 

u 

V 

V 

CO 


temperature,  K 
thrust,  N 

unit  veetor  tangential  to  a  surfaee 
divergent  half  angle  of  the  nozzle 
axial  veloeity,  m/s 
radial  veloeity  eomponent,  m/s 
veloeity  vector 
viscosity  exponent 


Subscripts 

0  =  plenum  or  stagnation 

c  =  curvature 

cont  =  continuum 

e  =  exit  condition 

fm  =  free  molecular 

0  =  orifice 

n  =  nozzle 

T]  =  value  normal  to  a  surface 

ref  =  reference  value 

t  =  throat  condition 

r  =  value  tangential  to  a  surface 

w  =  wall 


X  =  axial  component 
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Superscripts 

*  =  sonic  condition 

overbar=  mean  value 

prime  =  after  eollision  value 


Introduction 

Future  miero-  and  nanosatellites  will  require  low  thrust  propulsion  systems  for  orbital 
maneuvers  to  eounter  the  effects  of  trajectory  perturbations  and  for  fine  attitude  eontrol. 
Mieropropulsion  systems  for  these  vehieles  are  expeeted  to  have  thrust  in  the  range  from  1  pN  to 
1  mN.  For  propulsion  systems  that  utilize  nozzle  expansions,  the  low  thrust  levels  will  require 
either  low  pressure  operation  or  a  deerease  in  the  throat  dimensions,  both  potentially  leading  to 
low  Reynolds  number  flows.  The  further  development  of  effieient  miero-ehemioal  and 
eleetrothermal  thrusters  is  associated  with  high  stagnation  chamber  temperatures,  whieh  will  also 
eontinue  the  trend  toward  deereasing  the  Reynolds  number.  Therefore,  the  numerieal  and 
experimental  investigation  of  low  Reynolds  number  flows  has  become  inereasingly  important  for 
designing  effieient,  low  thrust  nozzles. 

It  is  well  known  that  the  Reynolds  number  is  a  measure  of  the  nozzle  effieiency  in  terms 
of  the  viscous  losses  inherent  in  the  subsonie  layer  near  nozzle  surfaees.  At  suffieiently  low 
Reynolds  numbers,  the  viscous  losses  within  a  mieronozzle  beeome  large  enough  for  the  speeifie 
impulse  along  the  nozzle  eenterline  to  deerease  from  the  nozzle  throat  to  the  exit  planed  perhaps 
making  the  divergent  nozzle  contour  useless.  In  the  transitional  flow  regime,  therefore,  a  thin- 
walled  sonie  orifiee  may  outperform,  or  at  least  perform  as  well  as,  a  typieal  mieronozzle. 
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Although  there  is  extensive  literature  related  to  rarefied  orifiee  flows,  there  have  been 
relatively  few  investigations  of  orifiee  performanee  in  terms  of  propulsive  parameters  sueh  as 
thrust  and  speeific  impulse.  ’  Investigations  of  low  Reynolds  number  flows  through  nozzles  of 
various  geometrie  shapes  are  numerous,  with  several  notable  eomputational  efforts  and 
experimental  studies^'^^  having  been  performed.  The  Navier-Stokes  (NS)  equations  with  veloeity 
slip  and  the  direet  simulation  Monte  Carlo  (DSMC)  method  are  typieally  used  to  model  low 
Reynolds  number  flows.  For  very  low  Reynolds  numbers  eorresponding  to  transitional  flow, 
reliable  results  ean  only  be  obtained  using  a  kinetic  approach.  A  previous  study  has  verified  that 
the  numerical  results  obtained  with  a  DSMC  flow  solver  reproduced  the  experimental  data  of 
Rothe^°  while  the  solution  of  the  NS  equations  with  no-slip  boundary  conditions  did  not.^ 
Numerical  simulations  referenced  above  have  shown  that  a  viscous  layer  occupies  a 
significant  fraction  of  a  nozzle’s  diverging  section  for  throat  Reynolds  numbers  less  than  about 
500.  As  the  Reynolds  number  further  decreases,  the  viscous  interactions  dominate  the  internal 
nozzle  flow  field,  limiting  nozzle  efficiency.  Because  an  orifice  with  a  small  thickness  to 
diameter  ratio  will  have  a  reduced  viscous  interaction  layer,  it  may  provide  an  adequate 
propulsive  efficiency  relative  to  a  nozzle  expansion  at  low  Reynolds  numbers. 

The  main  goal  of  this  work  is  to  study  experimentally  and  numerically  the  performance 
of  a  thin-walled,  squared  edge  orifice  and  a  DeLaval  nozzle  in  the  transitional  flow  regime.  The 
primary  objectives  are  to  compare  the  relative  performance  of  the  nozzle  and  orifice  geometries, 
find  the  range  of  Reynolds  numbers  where  the  gas  expansion  through  a  micronozzle  geometry  is 
no  longer  valid  from  an  efficiency  standpoint,  and  investigate  design  parameters  such  as 
expansion  angle  and  length  which  can  improve  micronozzle  performance  at  low  Reynolds 
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numbers.  The  study  also  extends  the  low  end  of  the  Reynolds  number  range,  eompared  to 
previous  studies,  by  extending  operating  eonditions  down  to  the  free  moleeule  flow  regime. 


Orifice  and  nozzle  thrust  for  an  inviscid  flow 

To  eompare  mieronozzle  and  orifiee  performanee,  as  preseribed  above,  it  is  advantageous 
to  use  the  relative  thrust  or  Isp  ratios  at  a  speeified  pressure  or  throat  Reynolds  number.  The 
Reynolds  number  at  a  nozzle  throat  or  a  thin-walled  orifiee  exit  plane  is  defined  as 

(1) 


The  throat  density,  p  ,  and  the  throat  veloeity,  u  ,  are  ealculated  from  measured  stagnation 
eonditions  and  the  quasi  one-dimensional  flow  relations.  The  viseosity  at  the  throat  or  orifiee  exit 
plane  is 


Thruster  performanee  is  determined  by  eonsidering  the  thrust  and  Isp  ratio  between  a 
nozzle  and  a  thin-walled  orifice  case  with  the  same  stagnation  conditions.  There  are  two  limits 
that  need  to  be  considered,  free  molecular  and  continuum.  The  free  molecule  analytical  solution 
for  an  infinitely  thin  orifice  case  is 

~^o,fm  ~  ■j/'oA?  (3) 
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There  is  no  analytieal  solution  available  for  the  free  molecule  nozzle  case,  since  it  depends 
inherently  on  the  nozzle  geometry  and  gas-surface  interactions. 

Assuming  there  is  no  background  pressure,  the  isentropic  continuum  thrust  from  an 
infinite  expansion  ratio  nozzle  is 


n,cont  ~  y\ 


f  ^  ll 

"  2  ' 

U-iJ 

[r  +  ij 

r+h 


A/ 


Po^t 


(4) 


Between  these  two  limits,  the  computational  thrust  is  calculated  for  zero  background 
pressure  as 

~  Pe^e  ~  Pe(^e  (^) 

This  expression  is  usually  reduced  to  3  =  (PeUg^  +  Pe)^e  finite  Reynolds  numbers  with 
supersonic  conditions,  which  cannot  be  used  for  the  nearly  sonic  cases  considered  here. 
Therefore,  Eq.  5  is  utilized  for  thrust  calculations  in  this  work.  To  calculate  the  thrust  from  a 
nozzle  or  orifice  from  Eq.  5  using  a  NS  solver,  one  needs  to  implement  the  formulation  of  the  NS 
equations  with  separate  temperature  components. 


Experimental  Setup 

All  thrust  measurements  were  performed  on  the  nano-Newton  Thrust  Stand  (nNTS) 
which  has  been  described  in  detail  by  Jamison,  et.  al.'"^  The  nNTS  was  installed  in  Chamber  IV 
of  the  Collaborative  High  Altitude  Elow  Eacilities  (CHAEE-IV)  which  is  a  3  m  diameter  by  6  m 
long  cylindrical,  high  vacuum  chamber.  The  facility  was  pumped  with  aim  diameter  diffusion 
pump  with  a  pumping  speed  of  42,000  L/s  for  helium  and  25,000  E/s  for  molecular  nitrogen. 
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The  ultimate  facility  pressure  was  approximately  10'^  Torr  with  all  operational  pressures  below 
10'"^  Torr.  A  previous  study has  shown  that  at  these  background  pressures  and  corresponding 
thrust  levels  there  is  no  evidence  of  background  pressure  effects  on  the  thrust  measurements  in 
CHAFF-IV. 


The  thin-walled  orifice  used  in  this  study  is  shown  schematically  in  Fig.  1.  The  1  mm 
diameter  orifice  was  attached  to  a  plenum  with  a  cross-sectional  area  much  larger  than  the  orifice 
area  to  help  ensure  uniform  flow.  The  orifice  was  machined  in  a  0.015  mm  thick  tantalum  shim 
giving  t/J  =  0.015.  The  thin-walled  geometry  was  chosen  to  minimize  the  viscous  effects  in  the 
orifice  flow. 

The  conical  De  Laval  nozzle  used  in  this  study  is  shown  schematically  in  Fig.  2.  The 
conical  nozzle  was  scaled  from  the  geometry  used  by  Rothe.'°  The  scaled  Rothe  geometry  has  a 


30°  subsonic  section,  a  relatively  sharp  1mm  diameter  throat  with  radius  of  curvature  a 


20°  diverging  section,  and  an  expansion  ratio  of  62.4.  This  geometry  was  selected  due  to  the 
extensive  experimental  data  that  exists'**  which  was  previously  used  to  verily  the  DSMC  model’s 
accuracy. '  The  nozzle  was  machined  from  aluminum  and  attached  to  a  cylindrical  aluminum 
plenum  and  mounted  on  the  nNTS.  Figure  3  shows  a  scanning  electron  microscope  image  of  the 
nozzle  side  wall  where  the  surface  features  caused  by  the  machining  process  are  clearly  evident. 
The  effect  of  the  rough  diverging  section  walls  on  the  nozzle’s  performance  will  be  investigated 
in  the  following  sections. 

The  stagnation  pressures  for  the  orifice  and  nozzle  were  measured  inside  their  respective 


plenums  using  calibrated  pressure  transducers.  The  propellant  was  introduced  to  the  plenum 
through  an  adjustable  needle  valve  located  downstream  of  a  mass  flow  meter.  In  the 
experimental  configuration,  the  mass  flow  meters  were  operated  in  the  continuum  regime 
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throughout  the  pressure  range  investigated.  The  propellants  used  were  molecular  nitrogen  and 
helium.  In  this  study,  the  stagnation  pressures  ranged  from  several  milliTorr  to  approximately  17 
Torr  for  both  propellants,  and  the  stagnation  temperature  was  measured  to  be  295  K.  The 
combination  of  stagnation  pressure  and  temperature  gave  maximum  Reynolds  numbers  of  350 
for  nitrogen  and  about  150  for  helium. 

The  nNTS  was  calibrated  using  the  gas  dynamic  and  electrostatic  calibration  techniques 
described  by  Selden  and  Ketsdever,'^  which  have  been  validated  through  experimental  and 
numerical  analysis.  A  unique  feature  of  the  nNTS  is  its  ability  to  measure  the  force  levels  of  the 
1.0  mm  orifice  and  nozzle  geometries  from  the  free  molecule  through  continuum  range.  The  low 
thrust  measuring  capability  of  the  nNTS  allowed  for  the  investigation  of  the  transitional  flow 
regime  overlooked  in  previous  low  Reynolds  number  nozzle  studies. 

Numerical  method 

The  flow  and  geometric  conditions  of  the  experimental  setup  described  above  were 
simulated  and  used  for  validation  of  the  DSMC  models  and  parameters.  All  walls  were  assumed 
to  be  at  a  temperature  of  295K  (as  measured  in  the  experiments),  and  the  gases  considered 
include  helium  and  nitrogen  at  a  stagnation  temperature  of  295K.  A  background  pressure  of  zero 
was  assumed  in  the  calculations.  The  throat  Reynolds  numbers  ranged  from  0.02  to  270  to 
adequately  cover  the  range  of  the  experimental  data.  After  the  numerical  model  was  analyzed,  it 
was  used  to  study  the  impact  of  nozzle  parameters  of  interest.  First,  the  half  angle  of  the 
diverging  section  of  the  nozzle  was  varied  from  the  experimental  case  {0=  20°)  to  25°,  30°,  and 
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40°  while  keeping  the  wall  length,  Iw,  eonstant.  Seeond,  the  nozzle  diverging  seetion  length  was 
then  shortened  from  the  experimental  length  to  L/2  and  L/3. 

The  SMILE  eomputational  tool  based  on  the  DSMC  method  was  used  in  the  numerieal 
part  of  this  work.  Details  on  the  tool  may  be  found  elsewhere.*’  The  SMILE  eapabilities  that 
were  used  in  the  present  work  inelude  two-level  reetangular  grids  adaptive  to  flow  gradients, 
different  grids  for  eollisions  and  maeroparameters,  and  parallel  implementation  with  a  statie  load 
balaneing  teehnique.  The  majorant  frequeney  seheme  was  employed  for  modeling  moleeular 
eollisions,  and  the  variable  hard  sphere  model  was  used  for  modeling  intermoleeular 
interaetions.  The  diserete  Larsen-Borgnakke  model*^’’**  with  temperature-dependent  rotational 
relaxation  number  was  utilized  for  rotation-translation  energy  transfer.  The  Maxwell  model  with 
different  values  of  aeeommodation  eoeffieients  was  used  to  model  gas-surfaee  interaetion.  The 
baekground  eell  size  was  a  fraetion  of  the  mean  free  path  at  the  throat  eonditions.  Over  2  million 
simulated  partieles  were  used  in  the  nozzle  eases,  and  over  6  million  partieles  were  used  in  the 
orifiee  eases. 

The  diffieulty  with  ealeulations  involving  subsonie  boundary  eonditions  is  attributed  to 
the  upstream  influenee  of  the  downstream  flow.  To  deal  with  these  effeets,  results  must  be 
independent  of  the  domain  size.  Eor  the  orifiee,  the  simulated  plenum  (subsonie  seetion)  was 
varied  in  the  axial  direetion  from  3  to  9  d,  and  in  the  radial  direetion  from  4  to  8  with  an 
exit  plane  Reynolds  number  of  200  to  show  that  the  maero-parameters  at  the  exit  are  independent 
of  the  plenum  size.  The  simulated  plenum  size  was  then  ehosen  to  be  9  x  8  for  all  oases.  The 
eomputational  domain  beyond  the  orifiee  exit  plane  (simulated  vaouum  faoility  domain)  was 
varied  in  the  axial  direetion  from  lid,  to  \5  d^  and  in  the  radial  direetion  from  5  d,  to  30  to 
show  the  independenoe  of  the  results  relative  to  the  vaouum  faoility  domain  size.  The  simulated 
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vacuum  facility  domain  was  then  chosen  to  be  W  d,  xSd^.  An  equilibrium  distribution  with 
values  of  stagnation  pressure  and  temperature  and  zero  average  velocity  was  used  at  the  top  and 
left  boundaries  of  the  plenum  in  all  orifice  calculations. 

For  the  nozzle,  the  computational  domain  beyond  the  exit  plane  (exhausting  into  the 
vacuum  facility)  was  varied  in  the  axial  direction  from  2  d,  to  20  d,  and  in  the  radial  direction 
from  \  d,  io  W  di  above  the  exit.  The  effect  of  having  a  larger  downstream  region  was  examined 
and  determined  to  have  negligible  effect.  This  exit  region  was  then  chosen  to  go  2  d,  beyond  the 
exit  and  1  d,  above  the  exit.  To  study  the  effect  of  the  subsonic  inflow  condition,  computations 
were  performed  for  a  plenum  of  5  x5d^  and  for  inflow  boundaries  loeated  at  M  =  0.1  and  0.3. 

The  effect  on  the  nozzle  thrust  with  a  throat  Reynolds  number  of  200  was  found  to  be  negligible. 

* 

The  computed  thrust  for  Re  =  60  was  found  to  be  up  to  3%  different  for  the  subsonic 
boundaries  loeated  at  M  =  0.1  and  0.3  compared  to  thrust  eomputed  with  the  plenum  domain 
included.  The  difference  was  due  to  the  impact  of  the  viscous  interaction  layer  in  the  eonverging 
section.  Therefore,  the  plenum  geometry  was  used  throughout  the  remainder  of  this  study,  and 
the  plenum  dimensions  remained  constant  for  all  subsequent  simulations. 

Results  and  Discussion 

A.  Orifice  and  nozzle  flow  structure 

First,  the  flow  structure  for  the  two  geometries  under  consideration  with  different  gas 
densities,  or  Reynolds  numbers,  is  examined  numerically.  Note  that  the  results  in  this  section  are 
shown  for  helium  flows  only,  since  the  general  flow  structure  for  nitrogen  is  qualitatively  similar 
for  the  same  Reynolds  number.  The  predicted  Mach  number  contours  for  the  orifiee  geometry 
with  Re*  =10  and  200  are  presented  in  Fig.  4.  Only  a  fraction  of  the  eomputational  domain  is 
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shown  here  in  order  to  illustrate  the  details  of  the  flow  elose  to  the  orifiee  plane.  In  the  figure,  the 
plenum  is  loeated  to  the  left  of  the  orifiee  plane,  and  the  vacuum  facility  domain  is  to  the  right. 

The  area  of  the  flow  impact  upstream  of  the  orifice  is  comparatively  small,  and  the  gas  is 
essentially  undisturbed  at  distances  larger  than  two  orifice  diameters.  The  impact  of  the  Reynolds 
number  on  the  subsonic  part  of  the  flow  (in  the  plenum)  is  visible  but  not  large,  with  the 
difference  in  the  Mach  number  amounting  to  about  40  percent  when  the  Reynolds  number 
increases  from  10  to  200.  The  flow  becomes  supersonic  at  a  small  distance  downstream  from  the 
orifice.  The  downstream  distance  of  the  sonic  line  increases  slightly  with  Reynolds  number.  A 
much  larger  number  of  collisions  in  the  plume  for  Re*  =  200  results  in  a  significantly  lower 
translational  temperature  downstream  from  the  orifice,  which  explains  the  higher  Mach  number 
values  when  compared  to  the  Re*  =10  case. 

The  impact  of  the  Reynolds  number  and  rarefaction  effects  is  more  pronounced  for  the 
nozzle  flow  because  of  gas-surface  interactions  resulting  in  fully  viscous  flow.  The  chamber  gas 
mean  free  path  for  Re*  =  10  is  about  10  ,  which  gives  a  Knudsen  number  of  approximately 

0.1  based  on  the  throat  diameter.  The  gas  rarefaction  for  this  Knudsen  number  is  high  enough  for 
the  nozzle  surface  to  influence  the  gas  in  the  diverging  and  converging  sections.  This  is  clearly 
seen  in  Fig.  5  where  the  Mach  number  contours  are  shown  for  Re*  =10  and  200.  The  surface 
effects  in  the  diverging  nozzle  section  result  in  subsonic  flow  in  a  majority  of  the  section,  with  a 
Mach  number  at  the  nozzle  exit  of  about  1.3. 

For  the  largest  Reynolds  number  considered.  Re*  =  200,  the  nozzle  flow  is  also 
dominated  by  surface  interactions.  In  this  case  however,  an  inviscid  core  starts  to  form,  and  the 
sonic  line  is  located  close  to  the  nozzle  throat.  Still,  the  viscous  interaction  layer  is  very  large  and 
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visible  even  in  the  converging  section.  At  the  nozzle  exit,  the  Mach  number  drops  from  about  3.2 
at  the  nozzle  exit  to  less  than  unity  at  the  lip. 

Relatively  large  mean  free  paths  and  strong  surface  effects  cause  a  significant 
nonequilibrium  between  translational  modes  in  both  the  orifice  and  the  nozzle  flows.  The  extent 
of  flow  nonequilibrium  is  illustrated  in  Fig.  6,  where  the  ratio  of  (the  temperature  calculated 
using  only  particle  velocities  in  the  axial  direction)  to  the  overall  translational  temperature  T  is 
presented  for  the  orifice  geometry,  ft  is  seen  that  for  Re*  =10  there  is  a  significant  degree  of 
translational  mode  nonequilibrium  inside  the  plenum.  In  the  core  flow  downstream  from  the 
orifice,  the  temperature  in  the  axial  direction  is  significantly  larger  than  that  in  the  radial 
direction  due  to  the  flow  expansion  (the  axial  component  of  velocity  is  largest  closer  to  the  axis). 
The  radial  component  is  larger  for  the  flow  angles  larger  than  about  45°  from  the  axis. 

For  a  nozzle  flow  expanding  into  a  vacuum,  a  strong  thermal  nonequilibrium  in  the 
region  of  the  nozzle  lip  and  in  the  back  flow  region  is  usually  expected.  For  rarefied  flows 
examined  in  this  work,  the  gas  is  in  thermal  nonequilibrium  even  inside  the  converging  section 
of  the  nozzle  as  shown  in  Fig.  7.  For  Re*  =10,  the  T^IT  ratio  inside  the  nozzle  reaches  a 
maximum  of  1.13  near  the  center  of  the  diverging  section  and  a  minimum  of  0.78  near  the  nozzle 
lip.  Thermal  nonequilibrium  is  less  significant  for  Re*  =  200.  Nevertheless,  the  minimum  value 
inside  the  nozzle  of  0.7  is  less  than  for  Re*  =10,  substantiating  the  use  of  a  kinetic  approach  for 
numerical  modeling  of  these  flows. 

A  more  quantitative  illustration  of  the  rarefaction  effect  of  the  flow  at  the  orifice  plane  is 
given  in  Fig.  8,  where  the  profiles  of  the  axial  velocity  are  presented  for  Re*  =  10,  60  and  200. 
Hereafter  7  =  0  corresponds  to  the  axisymmetric  centerline.  For  the  most  rarefied  flow. 

Re*  =10,  the  velocity  reaches  the  maximum  of  720  ml s  on  axis,  and  then  decreases  toward  the 
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lip;  however,  the  behavior  ehanges  qualitatively  for  larger  Reynolds  numbers.  For  Re*  =  60  and 
Re*  =  200,  there  is  a  maximum  at  about  Y  =  0.4  mm,  and  the  veloeity  at  the  axis  is  10  to  20 
pereent  lower  than  the  maximum  value.  The  maximum  in  the  axial  veloeity  is  due  to  the  low 
density  expansion  region  above  the  plenum  (large  angles  from  the  eenterline).  The  moleeules 
from  that  region  have  a  very  low  probability  of  baekseattering  to  the  plenum  eompared  to  the 
moleeules  that  leave  the  orifiee  near  the  eenterline.  Therefore,  the  impaet  of  the  returned 
moleeules  is  lower  at  larger  radial  distanees. 

The  axial  veloeity  profiles  aeross  the  nozzle  exit  are  shown  in  Fig.  9.  It  is  interesting  to 
note  that  whereas  the  veloeity  is  mueh  larger  on  the  nozzle  axis  for  Re*  =  200  than  for  Re*  =10, 
the  veloeity  is  only  weakly  dependent  on  the  Reynolds  number  at  the  nozzle  lip.  For  all  oases 
under  oonsideration,  the  axial  veloeity  tends  to  deorease  with  radial  distanoe. 

B.  Comparison  of  experimental  and  numerioal  results 

The  oomparative  performanoe  analysis  of  an  orifiee  and  a  nozzle  has  been  oonduoted 
numerioally  and  experimentally  at  different  oonditions  in  the  free  moleoular  and  transitional  flow 
regimes.  The  oomputational  results  for  helium  flow  over  the  entire  range  of  Reynolds  numbers 
oonsidered  are  given  in  Table  1.  As  expeoted,  thrust  and  speoifio  impulse  inorease  with 
inoreasing  Reynolds  number  for  both  the  orifiee  and  nozzle.  The  rate  of  this  inorease,  however,  is 
different  for  these  geometries.  As  a  result,  the  relative  performanoe,  represented  by  the  ratio  of 
the  nozzle  to  orifiee  thrusts  and  speoifio  impulses,  has  a  minimum. 

Consider  the  meohanisms  that  govern  the  relative  performanoe  oharaoteristios  for 
different  Reynolds  numbers.  At  large  Reynolds  numbers,  where  the  flow  is  near-oontinuum  and 
the  surfaoe  effeots  are  less  pronounoed,  the  thrust  ratio  should  approaoh  some  oontinuum  limit 
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value.  The  increase  in  flow  rarefaction  causes  a  larger  impact  of  the  surface,  and  the  surface 
effect  is  more  important  in  the  nozzle  geometry,  where  the  surface  to  volume  ratio  is  much  larger 
than  for  the  orifice. 

In  the  opposite  limit  of  the  free  molecular  flow,  the  nozzle  geometry  and  the  gas-surface 
interaction  law  primarily  influence  the  thrust  ratio.  It  may  generally  be  much  larger  than  unity  for 
a  specular  surface,  and  less  than  unity  for  a  diffuse  surface.  In  the  presence  of  molecular 
collisions,  when  the  flow  becomes  transitional,  the  thrust  ratio  is  expected  to  decrease  as  the 
viscous  interaction  layer  inside  the  nozzle  becomes  important.  As  a  result,  a  minimum  is 
observed  between  the  free  molecular  and  continuum  limits.  An  interesting  fact  is  that  the 
minimum  is  located  about  a  Reynolds  number  of  10  for  the  thrust  ratio,  and  60  for  the  specific 
impulse  ratio. 

The  orifice  and  nozzle  thrust  in  the  free  molecular  and  early  transitional  regime  has  been 
examined  in  more  detail.  The  obtained  numerical  results  have  been  compared  with  the 
experimental  measurements  in  Figures  10  and  11,  which  present  a  comparison  for  a  helium  flow 
inside  the  orifice  and  the  nozzle,  respectively,  in  the  range  of  Reynolds  numbers  from  0.02  (free 
molecular  regime)  to  23.  Note  that  the  orifice  results,  both  experimental  and  computed,  begin  to 
deviate  from  free  molecule  theory  at  Reynolds  numbers  greater  than  one  and  remain  lower  than 
the  corresponding  sonic  nozzle  continuum  value  at  Re  >  20. 

The  experimental  and  numerical  results  show  that  the  nozzle  thrust  does  not  yet  approach 
its  continuum  value  by  a  Reynolds  number  of  20,  with  the  difference  of  about  a  factor  of  two  as 
shown  in  Fig.  11.  The  difference  between  the  measurements  and  the  calculations  is  somewhat 
larger  for  the  nozzle  than  for  the  orifice,  with  the  experimental  slope  being  steeper.  This 
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difference  is  attributed  primarily  to  the  gas-surface  interaction  model  used  in  the  computations  as 
discussed  in  more  detail  below. 

Comparison  of  the  numerical  and  experimental  results  in  a  wider  range  of  Reynolds 
numbers  is  presented  in  Fig.  12,  where  the  nozzle-to-orifice  thrust  ratio  is  given  for  a  helium 
flow.  Comparison  of  the  experimental  data  with  the  modeling  results  for  fully  diffuse  walls 
shows  that  the  calculations  are  in  good  agreement  with  the  data  at  larger  Reynolds  numbers  and 
overpredict  it  at  lower  pressures.  The  numerical  error  due  to  statistical  scatter,  convergence 
issues,  and  boundary  conditions  is  estimated  to  be  on  the  order  of  one  percent  in  a  single  thrust 
value  and  two  percent  in  the  thrust  ratio.  The  reason  for  the  difference  in  the  slope,  apart  from 
experimental  uncertainties,  lies  therefore  in  the  physical  model  used  in  the  DSMC  calculations. 
The  largest  uncertainty  in  the  physical  model  is  related  to  the  gas-surface  interaction  model.  It 
also  has  the  biggest  impact  on  the  thrust  slope  compared  to  the  molecular  interaction  potential, 
chamber  conditions,  or  facility  background  gas  pressure.  A  number  of  computations  have 
therefore  been  performed  for  incomplete  energy  and  momentum  accommodation. 

It  is  suggested  by  Arkilic,  et  al.  that  a  tangential  momentum  accommodation  coefficient 
of  0.8  fits  best  for  low  Reynolds  number  flows.  The  results  for  the  Maxwell  model  with  80 
percent  diffuse  and  20  percent  specular  reflection  are  given  in  Fig.  12.  Note  that  the  energy 
accommodation  coefficient  is  not  very  important  since  the  wall  temperature  is  equal  to  the 
stagnation  temperature.  An  addition  of  surface  specularity  results  in  larger  thrust  for  the  nozzle 
flow  in  the  entire  range  of  Reynolds  numbers.  The  general  trend  of  the  numerical  points  is  even 
flatter  than  for  the  fully  diffuse  surface  accommodation,  especially  for  low  Reynolds  numbers 
where  it  strongly  overpredicts  the  experimental  values. 
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A  close  examination  of  the  surfaee  strueture  inside  the  aetual  nozzle  manifested  a  very 
rough,  groove-like  strueture,  as  illustrated  in  Fig.  3,  with  mieron-size  grooves  set  out 
perpendieular  to  the  main  flow  direetion.  The  evident  surfaee  roughness  prompted  the  authors  to 
introduee  a  simple  model  to  simulate  the  roughness.  In  this  model,  a  fraetion  of  partieles,  a,  is 
refleeted  diffusely  on  the  surfaee,  and  a  eomplementary  fraetion  of  partieles,  (1  -  a),  is  refleeted 
with  veloeities  both  tangential  and  normal  to  the  surfaee,  and  respeetively,  taking  an 

opposite  sign  as  eompared  to  the  ineident  partiele  veloeities.  This  type  of  refleotion  is  hereafter 
referred  to  as  antispeeular.  The  diffuse-antispeeular  refleetion  kernel  may  be  written  as 
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The  results  of  eomputations  with  «  =  80  (80  pereent  diffuse  and  20  pereent  antispeeular) 
are  also  shown  in  Fig.  12.  It  is  elearly  seen  that  whereas  the  points  generally  fall  further  from  the 
experimental  data  than  for  the  fully  diffuse  model,  the  diffuse-antispeeular  model  prediets  a 
thrust  ratio  slope  elose  to  the  experimental  one. 

The  eomparison  of  the  experimental  and  eomputed  nozzle-to-orifiee  thrust  ratios  for 
moleeular  nitrogen  is  presented  in  Fig.  13.  The  general  trend  observed  for  the  helium  flow, 
where  the  numerieal  results  were  elose  to  the  measurements  for  the  moderate  Reynolds  numbers, 
holds  for  the  nitrogen  flow  as  well. 

C.  Impaet  of  the  nozzle  geometry 
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Previous  numerical  studies^  of  low  Reynolds  number  nozzle  flows  analyzed  the  variation 
of  the  propulsive  efficiency  for  different  nozzle  contours:  conical,  bell  and  trumpet.  Zelesnik,  et 
al^  showed  that  the  trumpet-shaped  nozzle  has  about  5%  higher  efficiency  compared  to  a  conical 
nozzle  for  a  stagnation  temperature  of  approximately  300K.  However,  the  easier  to  manufacture 
conical  nozzle  was  shown  to  have  a  higher  mass  discharge  coefficient. 

In  an  attempt  to  define  an  optimum  conical  diverging  nozzle  configuration  in  the  low 
Reynolds  number  regime,  a  set  of  computations  has  been  performed  for  different  geometrical 
configurations,  with  varied  nozzle  length  and  diverging  section  angle.  The  results  were  obtained 
for  helium  test  gas  and  Reynolds  numbers  of  60  and  200.  The  surface  accommodation  in  these 
computations  was  assumed  to  be  fully  diffuse. 

To  better  understand  the  impact  of  the  viscous  effects  on  the  nozzle  efficiency  for 
different  nozzle  geometries  and  chamber  pressures,  the  specific  impulse  has  been  calculated  as  a 
function  of  the  axial  distance  from  the  nozzle  throat  to  the  exit.  Figure  14a  shows  the  specific 
impulse  profiles  for  three  values  of  the  nozzle  length  for  Re*  =  60 .  Consider  first  the  baseline 
length  of  L  =  1.07  x  lO^^m .  Figure  14a  shows  that  the  specific  impulse  increases  sharply  over  the 
first  millimeter  past  the  throat.  Further  downstream  the  specific  impulse  had  a  maximum  and 
then  decreased  toward  the  nozzle  exit  plane  as  the  influence  of  the  viscous  layer  increases.  Note 
that  the  maximum  that  occurs  at  X  «  0.4T  was  about  5%  larger  than  that  at  the  nozzle  exit. 

The  numerical  results  indicated  that  shortening  the  nozzle  generally  resulted  in  higher  axial  flow 
velocities  at  the  nozzle  exit  plane  close  to  the  centerline  due  to  a  smaller  impact  of  the  viscous 
layer.  However,  the  total  1^^  increase  amounts  to  only  about  4%  when  the  nozzle  was  three 

times  shorter  than  the  baseline  case  for  Re*  =  60 . 
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As  expected,  an  increase  in  the  Reynolds  number  resulted  in  better  nozzle  efficiencies,  as 
shown  in  Fig.  14b  for  Re*  =  200  and  three  values  of  the  nozzle  length.  Still,  there  was  a 
maximum  in  the  specific  impulse  observed  inside  the  nozzle.  The  location  of  the  maximum 
shifted  downstream  as  compared  to  Re*  =  60  due  to  the  thinner  viscous  layer.  The  value  of  at 

the  nozzle  exit  was  about  2%  lower  than  the  maximum  value  inside  the  nozzle.  The  shortening  of 
the  nozzle  causes  about  a  1%  decrease  with  the  relative  position  of  the  maximum  being  closer  to 
the  exit  plane  for  shorter  nozzles.  A  smaller  influence  of  the  viscous  layer  for  larger  Reynolds 
number  flows  was  the  reason  for  an  insignificant  impact  of  the  nozzle  length.  The  difference 
between  the  cases  of  L/2  and  L/3  in  terms  of  the  nozzle  efficiency  was  less  than  1%. 

The  second  nozzle  geometry  parameter  varied  was  the  angle  of  the  diverging  section,  ft  is 
well  known  that  conventional  high  Reynolds  number  conical  nozzles  are  usually  designed  with  a 
half-angle  of  about  15°,  where  the  nozzle  efficiency  is  maximum.  For  low  Reynolds  numbers,  the 
viscous  effects  become  a  significant  issue  that  degrades  nozzle  performance.  Therefore,  it 
appeared  reasonable  to  increase  the  half-angle  of  the  diverging  section  in  order  to  improve  the 
efficiency.  The  numerical  studies  in  this  work  were  carried  out  for  the  half-angle  from  20° 
(baseline  geometry)  up  to  40°.  The  angle  was  increased  while  maintaining  a  constant  diverging 
section  wall  length,  Iw. 

Specific  impulse  as  a  function  of  axial  station  for  four  different  half-angle  values  and 
Re*  =  60  is  given  in  Fig.  15a.  While  the  general  shape  of  the  1^^  distribution  does  not  change 

significantly  with  angle,  there  were  several  important  effects  associated  with  the  angle  increase. 
In  contrast  to  the  nozzle  length  change  at  a  constant  angle,  where  the  flow  was  essentially 
constant  over  the  first  millimeter  from  the  throat,  the  varying  angle  impacted  the  specific  impulse 
immediately  after  the  throat.  A  larger  angle  caused  higher  values  of  axial  velocity,  which 
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generally  resulted  in  higher  thrust  and  Isp  in  the  diverging  seetion  of  the  nozzle.  Further 
downstream,  the  viseous  layer  deereased  the  nozzle  effieieney,  whieh  was  espeeially  pronouneed 
for  smaller  angles.  A  notieeably  better  nozzle  performanee  was  seen  for  larger  angles.  Every  5° 
inerease  in  the  angle  resulted  in  about  a  2%  inerease  in  the  speeifie  impulse  at  the  nozzle  exit 
plane  for  Re*  =  60 .  At  larger  Reynolds  numbers,  the  inerease  of  the  nozzle  angle  does  not  have  a 
strong  impaet  on  the  flow  in  the  region  of  the  throat,  as  illustrated  in  Fig.  15b. 

The  thrust  and  speeifie  impulse  results  for  all  geometrie  eonditions  are  summarized  in 
Table  2.  The  general  eonelusion  that  ean  be  drawn  from  Table  2  is  that  the  inerease  of  the  nozzle 
angle  and/or  deerease  of  the  nozzle  length  inereased  the  thrust  and  improved  the  nozzle 
effieieney  for  Reynolds  numbers  less  than  approximately  100.  This  geometrie  improvement 
deereased  with  inereasing  Reynolds  number,  and  is  within  one  or  two  pereent  for  all  geometries 
eonsidered  at  Re*  =  200 . 


Numerical  and  Experimental  Uncertainties 

There  are  several  uneertainties  assoeiated  with  the  experimental  effort  in  this  study.  First, 
eorreetions  to  the  foree  measurements  due  to  the  faeility  baekground  gas  have  not  been 
ineorporated  into  the  data  analysis.  However,  a  previous  study^^  indieated  that  the  thrust  may  be 
effeeted  by  less  than  0.5%  at  the  experimental  eonditions  of  this  work.  Seeond,  the  standared 
deviation  of  the  stagnation  pressure  and  mass  flow  measurements  was  found  to  be  within 
approximately  1%  over  the  range  used.  Third,  thrust  stand  ealibration  of  defieetion  versus 
applied  foree  has  been  approximated  to  be  within  3%.  For  a  given  applied  foree  to  the  thrust 
stand,  the  standard  deviation  of  the  stand’s  defieetion  was  less  than  1%;  however,  the  aeeuraey 
of  the  ealibration  system  must  also  be  taken  into  aeeount.  Finally,  there  was  some  error 
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associated  with  the  manufacturing  of  the  orifice  and  nozzle.  In  general,  the  nozzle  throat 
diameter  is  known  to  within  0.5%.  Because  the  orifiee  was  machined  from  a  thin  tantalum  shim, 
the  measurement  of  its  diameter  is  not  quite  as  accurate  and  is  estimated  to  be  within 
approximately  2%. 

In  addition  to  the  experimental  uncertainties,  there  are  a  number  of  numerical 
uncertainties.  Grid  resolution,  the  maximum  number  of  simulated  molecules,  effects  of  the 
subsonic  boundary  conditions,  and  the  gas-gas  collision  models  all  account  for  a  numerical 
uncertainty  estimated  to  be  on  the  order  of  1  to  2%.  The  effects  of  the  nozzle  surface  roughness 
present  in  the  experimental  device  are  not  known,  although  a  simple  diffuse-anti-specular  model 
was  used.  However,  it  is  not  known  to  what  extent  the  nozzle  surface  roughness  contributes  to 
uncertainties  in  either  the  numerical  or  experimental  data. 


Conclusions 


Low  Reynolds  number  flows  through  nozzle  and  orifice  geometries  have  been 
investigated  both  numerically  and  experimentally.  The  DSMC  results  show  that  the  viscous 
effects  dominate  the  flow  for  Re*  =  60.  This  result  has  been  seen  in  previous  studies^  and 
prompted  the  investigation  of  thin-walled  orifices  for  thrust  generation,  since  presumably  an 
orifice  would  not  exhibit  large  viscous  interaction  regions.  Experimental  results  indicate  that  at 
Re*  <100  the  thrust  generated  by  an  orifice  is  indeed  larger  than  that  of  the  nominal  conical 
nozzle  for  helium  and  nitrogen.  However,  the  nozzle  remains  more  efficient  (based  on  fip)  than 
the  orifice  throughout  the  range  of  kp  investigated.  There  is  a  minimum  in  the  nozzle  to  orifice 
specific  impulse  ratio  around  Re*  =  60  for  helium.  Therefore  at  some  Reynolds  numbers,  it  may 
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be  advantageous  to  use  a  simple  orifiee  geometry  in  this  range  over  a  nozzle  geometry  that  is 
eomplieated  to  fabrieate. 

In  general,  the  experimental  and  numerieal  results  for  the  nozzle  to  orifiee  thrust  ratio  are 
in  good  agreement  (within  a  few  pereent).  The  results  obtained  with  an  anti-speeular  gas-surfaee 
interaeton  model  for  DSMC  show  better  qualitative  agreement  with  the  slope  of  the  thrust  ratio 
data  than  the  standard  Maxwell  model.  The  anti-speeular  model  was  investigated  based  on 
seanning  eleetron  mieroseope  images  of  the  nozzle’s  surfaee  roughness  eaused  by  the  maehining 
proeess.  Although  the  exaet  influenee  of  the  surfaee  roughness  is  unknown,  it  appears  to  have  a 
measurable  effeet  on  the  nozzle’s  performanee. 

Onee  the  DSMC  model  results  were  verified  using  experimental  data,  the  model  was  used 
to  assess  geometrieal  variations  on  the  eonieal  nozzle’s  performanee.  Both  the  diverging  seetion 
length  and  angle  were  varied  to  assess  the  effeet  on  nozzle  performanee.  Little  effeet  was  seen  in 
improving  the  performanee  for  Re*  =  200  where  all  variations  were  within  2%.  However,  an  8% 
improvement  in  fip  was  seen  for  a  Re*  =  60  helium  flow  when  the  expansion  angle  was  doubled 
from  20°  to  40°.  At  Re*  =  60,  it  was  also  shown  that  shortening  the  diverging  seetion  length 
improved  the  nozzle’s  performanee.  Previous  studies^  have  shown  that  using  trumpet  geometries 
may  further  improve  the  performanee  over  a  eonieal  geometry  by  as  mueh  as  5%  at  low 
Reynolds  number.  Therefore,  reasonable  improvement  of  nozzle  performanee  for  Re*  <100  ean 
be  realized  through  the  use  of  numerieal  simulations  to  optimize  geometry. 

Acknowledgments 


22 


This  work  was  supported  in  part  by  the  U.S.  Air  Force  Office  of  Scientific  Research  and 
the  Propulsion  Directorate  of  the  Air  Force  Research  Laboratory  at  Edwards  Air  Force  Base, 
California.  The  authors  wish  to  thank  Mr.  Andrew  Jamison  and  Mr.  Taylor  Lilly  for  their 
invaluable  assistance  with  this  research. 


References 

'ivanov,  M.,  Markelov,  G.,  Ketsdever,  A.,  and  Wadsworth,  D.,  “Numerical  Study  of  Cold 
Gas  Micronozzle  Flows,”  AIAA  Paper  99-0166,  January  1996. 

Ketsdever,  A.,  “Thrust  Measurements  of  an  Underexpanded  Orifice  in  the  Transitional 
Regime,”  Proceedings  of  the  23rd  International  Symposium  on  Rarefied  Gas  Dynamics,  edited 
by  A.  Ketsdever  and  E.P.  Muntz,  American  Institute  of  Physics,  New  York,  2002,  pp.  1057- 
1064. 

Alexeenko,  A.,  Gimelshein,  S.,  Levin,  D.,  Ketsdever,  A.,  and  Ivanov,  M., 

“Measurements  and  Simulation  of  Orifice  Elow  for  Micropropulsion  Testing,”  Journal  of 
Propulsion  and  Power,  Vol.  19,  No.  4,  2003,  pp.  588-594. 

"^Kim,  S.,  “Calculations  of  Low-Reynolds-Number  Resistojet  Nozzles,”  Journal  of 
Spacecraft  and  Rockets,  Vol.  31,  No.  2,  1994,  pp.  259-264. 

^Boyd,  L,  Jafry,  Y.,  and  Vanden  Beukel,  J.,  “Particle  Simulations  of  Helium 
MicxoihmsiQxVXosNS,"  Journal  of  Spacecraft  and  Rockets,  Yo\.  31,  No.  2,  1994,  pp.  211-211 . 

^Zelesntk,  D.,  Micci,  M.,  and  Long,  L.,  “Direct  Simulation  Monte  Carlo  Model  of  Low 
Reynolds  Number  Nozzle  Llows,”  Journal  of  Propulsion  and  Power,  Vol.  10,  No.  4,  1994,  pp. 
546-553. 


23 


’Chung,  C.-H.,  Kim,  S.,  Stubbs,  R.,  and  De  Witt,  K.,  “Low-Density  Nozzle  Flow  by  the 
Direet  Simulation  Monte  Carlo  and  Continuum  Mehtods,”  Journal  of  Propulsion  and  Power, 
Vol.  11,  No.  1,  1995,  pp.  64-78. 

g 

Alexeenko,  A.,  Levin,  D.,  Gimelshein,  S.,  Collins,  R.,  and  Reed,  B.,  “Numerieal  Study 
of  Flow  Strueture  and  Thrust  Performanee  for  3-D  MEMS-Based  Nozzles,”  AIAA  Paper  2002- 
3194,  June  2002. 

^Mureh,  C.,  Broadwell,  J.,  Silver,  A.,  and  Mareisz,  T.,  “Low- Thrust  Nozzle 
Performance,”  AIAA  Paper  68-091,  January  1968. 

'°Rothe,  D.,  “Eleetron-Beam  Studies  of  Viseous  Elow  in  Supersonie  lAozzles,”  AIAA 
Journal,  Vol.  9,  No.  5,  1971,  pp.  804-810. 

'' Whalen,  M.,  “Eow  Reynolds  Number  Nozzle  Elow  Study,”  NASA  Technical 

Memorandum  100130,  July  1987. 

12 

Eempert,  W.,  Boehm,  M.,  Jiang,  N.,  Gimelshein,  S.,  and  Eevin,  D.,  “Comparison  of 
Moleeular  Tagging  Veloeimetry  Data  and  Direet  Simulation  Monte  Carlo  Simulation  in 
Supersonie  Miero  Jet  Elows,”  Expt.  In  Fluids,  Vol.  34,  2003,  pp.  403-411. 

1 3 

Jamision,  A.,  and  Ketsdever,  A.,  “Eow  Reynolds  Number  Performanee  Comparison  of 
an  Underexpanded  Orifiee  and  a  DeEaval  Nozzle,”  Proceedings  of  the  23rd  International 
Symposium  on  Rarefied  Gas  Dynamics,  edited  by  A.  Ketsdever  and  E.P.  Muntz,  Ameriean 
Institute  of  Physies,  New  York,  2002,  pp.  557-564. 

'"^Jamison,  A.,  Ketsdever,  A.,  Muntz,  E.P.,  “Gas  Dynamie  Calibration  of  a  Nano-Newton 
Thrust  Stand,”  Review  of  Scientific  Instruments,  Vol.  73,  No.  10,  2002,  pp.  3629-3637. 


24 


'^Ketsdever,  A.,  "Facility  Effects  on  Performance  Measurements  of  Mieropropulsion 
Systems  Whieh  Utilize  Gas  Expansion,"  Journal  of  Propulsion  and  Power,  Vol.  18,  No.  4,  2002, 
pp.  797-804. 

'^Selden,  N.,  Ketsdever,  A.,  “Comparison  of  Eoree  Balanee  Calibration  Teehniques  for 
the  nano-Newton  Range,”  Review  of  Seientifie  Instruments,  Vol.  74,  No.  12,  2003,  pp.  5249- 
5254. 

'^Ivanov,  M.S.,  Markelov,  G.N.,  Gimelshein,  S.E.  “Statistieal  simulation  of  reaetive 

rarefied  flows:  numerieal  approaeh  and  applieations,”  AIAA  Paper  98-2669,  June  1998. 

1 8 

Ivanov,  M.S.,  Rogasinsky,  S.V.  “Analysis  of  the  numerieal  teehniques  of  the  direet 
simulation  Monte  Carlo  method  in  the  rarefied  gas  dynamies,”  Soviet  J.  Numer.  Anal.  Math. 
Modeling,  Vol.  3,  No.  6,  1988,  pp.  453-465. 

'^Borgnakke,  C.  and  Earsen,  P.  S.,  “Statistieal  eollision  model  for  Monte  Carlo 

simulation  ofpolyatomie  gas  mixture,”/.  Comp.  Phys.  18:405-420,  1975. 

20 

Gimelshein,  S.E.,  Boyd,  I.D.,  Ivanov,  M.S.  “Modeling  of  internal  energy  transfer  in 

plume  flows  of  polyatomie  moleeules  by  the  DSMC  method,”  AIAA  Paper  99-0738.  Reno,  NV, 

January  1999. 

2 1 

Arkilie,  E.  B.,  Breuer,  K.  S.,  and  Sehmidt,  M.  A.,  “Mass  flow  and  tangential  momentum 
aeeommodation  in  silieon  mieromaehined  ehannels,”  Journal  of  Fluid  Mechanics,  Vol.  437, 
2001,  pp.  29-43. 


25 


Table  and  Figure  Captions: 

Table  1:  Numerical  data  for  the  performance  of  the  nominal  orifice  and  nozzle  geometries  as  a 
function  of  Reynolds  number. 

Table  2:  Numerical  data  for  the  performance  of  various  nozzle  geometries. 

Figure  1 ;  Thin-walled  orifice  geometry. 

Figure  2:  Experimental  geometry  of  the  conical  DeLaval  nozzle. 

Figure  3;  SEM  imaging  showing  the  surface  roughness  of  the  expanding  section  of  the  conical 
nozzle. 

Eigure  4;  Predicted  Mach  number  contours  for  helium  orifice  flow  at  Re*  =10  (top)  and  200 
(bottom). 

Eigure  5;  Predicted  Mach  number  contours  for  helium  flow  inside  the  nozzle  at  Re*  =10  (top) 
and  200  (bottom). 

Eigure  6;  Predicted  Tx/T  contours  for  helium  orifice  flow  at  Re*  =10  (top)  and  200  (bottom). 
Eigure  7;  Predicted  Tx/T  contours  for  helium  flow  inside  the  nozzle  at  Re*  =10  (top)  and  200 
(bottom). 

Eigure  8;  Predicted  axial  velocity  profiles  along  the  orifice  plane  for  different  Reynolds  numbers 
for  a  helium  flow.  Y=0  corresponds  to  the  symmetry  axis. 

Eigure  9;  Predicted  axial  velocity  profiles  across  the  nozzle  exit  plane  for  different  Reynolds 
numbers  for  a  helium  flow.  Y=0  corresponds  to  the  symmetry  axis. 

Eigure  10;  Orifice  thrust  comparison  between  experimental  and  numerical  data  for  helium. 
Eigure  11;  Nozzle  thrust  comparison  between  experimental  and  numerical  data  for  helium. 
Eigure  12;  Ratio  of  the  nozzle  to  orifice  thrust  for  helium  as  a  function  of  Reynolds  number. 


26 


Figure  13;  Ratio  of  the  nozzle  to  orifice  thrust  for  nitrogen  as  a  function  of  Reynolds  number. 
Figure  14:  Impact  of  the  nozzle  length  on  performance  for  helium,  a)  Re*  =  60,  b)  Re*  =  200. 
Figure  15;  Impact  of  the  nozzle  expansion  angle  on  performance  for  helium,  a)  Re*  =  60,  b)  Re* 
=  200. 


27 


Re* 

Kn 

Po  [Pa] 

To  [N] 

Tn  [N] 

Tn/To 

ISPo  [S] 

ISPn  [S] 

ISPn  /  ISPo 

0.02 

48.58 

3.59E-01 

1.41E-07 

1.28E-07 

0.907 

102.48 

127.09 

1.24 

0.2 

4.883 

3.57E+00 

1.43E-06 

1.26E-06 

0.881 

102.91 

126.17 

1.23 

2 

0.488 

3.57E+01 

1.64E-05 

1.38E-05 

0.805 

108.23 

127.57 

1.18 

5 

0.195 

8.92E+01 

4.81E-05 

3.83E-05 

0.796 

114.19 

128.67 

1.13 

10 

0.098 

1.78E+02 

1.13E-04 

8.89E-05 

0.786 

122.78 

129.85 

1.06 

23 

0.042 

4.10E+02 

3.02E-04 

2.49E-04 

0.825 

128.96 

132.79 

1.03 

60 

0.016 

1.07E+03 

8.69E-04 

8.01E-04 

0.922 

138.42 

140.53 

1.02 

90 

0.011 

1.61E+03 

1.34E-03 

1.30E-03 

0.970 

141.34 

145.01 

1.03 

120 

0.008 

2.14E+03 

1.79E-03 

1.81E-03 

1.01 

143.61 

148.86 

1.04 

200 

0.005 

3.57E+03 

3.04E-03 

3.28E-03 

1.08 

144.62 

155.38 

1.07 

Table  1 
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Angle 

(deg) 

* 

Re 

Length 

Tn  (N) 

Ispn  (s) 

ISPn  /  ISPNom 

Nominal  L 

8.01E-04 

140.53 

1.00 

60 

lw/2 

8.25E-04 

144.14 

1.03 

20 

lw/3 

8.33E-04 

145.55 

1.04 

Nominal  L 

3.28E-03 

155.38 

1.00 

200 

lw/2 

3.29E-03 

157.14 

1.01 

lw/3 

3.29E-03 

157.23 

1.01 

60 

Iw 

8.33E-04 

144.14 

1.03 

25 

lw/2 

8.42E-04 

146.43 

1.04 

200 

Iw 

3.32E-03 

157.57 

1.01 

lw/2 

3.34E-03 

158.17 

1.02 

30 

60 

Iw 

8.58E-04 

148.14 

1.05 

200 

Iw 

3.33E-03 

159.12 

1.02 

40 

60 

Iw 

8.82E-04 

151.76 

1.08 

200 

Iw 

3.34E-03 

159.10 

1.02 

Table  2 
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Aluminum  Plenum 
65  mm  H  x  35  mm  W 


Figure  1 
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Figure  2 
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